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We investigate the stability of a family of spherically sym¬ 
metric static solutions in vacuum Brans-Dicke theory (with 
w = 0) recently described by van Putten. Using linear pertnr- 
bation theory, we find one exponentially growing mode for ev¬ 
ery member of the family of solutions, and thus conclude that 
the solutions are not stable. Using a previously constructed 
code for spherically symmetric Brans-Dicke, additional evi¬ 
dence for instability is provided by directly evolving the static 
solutions with perturbations. The full non-linear evolutions 
also suggest that the solutions are black-hole-threshold criti¬ 
cal solutions. 
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Recently van Putten has proposed a one parameter 
family (parameter e) of solutions to spherically sym¬ 
metric Brans-Dicke theory for use in numerical rela¬ 
tivity as an approximate black hole [||. These solu¬ 
tions have the attractive property that for small val¬ 
ues of this parameter, the “exterior” solution approaches 
that of Schwarzschild. However, the event horizon of 
Schwarzschild is replaced with a high red-shift horizon 
and all metric components remain finite at this horizon. 
In addition to regularizing the horizon, these solutions 
have a global timelike coordinate. 

As van Putten stresses, these solutions could have 
promise for numerical relativity because of the difficul¬ 
ties that arise when dealing numerically with boundary 
conditions at the horizon of a black hole. To be useful 
as approximate black holes, however, the solutions, like 
Schwarzschild, must be stable. 

In this Communication, we recompute the solutions 
considered by van Putten and carry out a linear pertur¬ 
bation analysis about them. In so doing, we find, for 
generic e, modes which grow exponentially in time. We 
also directly evolve perturbations on the background of 
these solutions and confirm the instability predicted by 
linear theory. Thus although of some theoretical inter¬ 
est, these solutions are unlikely to be of direct use in the 
context of mocking-up a black hole in general relativity. 

To begin, let us review the static family of solutions 
considered by van Putten. We note that they were first 
written down by Brans and Dicke ^ but in an isotropic 
coordinate system as opposed to the Schwarzschild-like 
coordinates that van Putten uses. 

We work in Brans-Dicke theory and assume spherical 
symmetry. We choose a coordinate system such that the 
metric has the form 

ds^ = + r^dUl^. ( 1 ) 

Using (fir, t) for the Brans-Dicke field, the field equations 
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where, in vacuum, the Brans-Dicke stress tensor is given 
by 

{<!)■,pu - guu^4 >), (3) 

and uj is the Brans-Dicke coupling constant . The field 
(j) satisfies the covariant wave equation 

n<j) = ^matter ^ (^4^ 
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This equation, along with van Putten’s restriction to 
w = 0, simplifies Eq. (^ to = (jj-fiv/Sn. As the 
stress tensor is traceless, the field equations may then be 
written in the simple form 
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where is the usual Ricci curvature tensor. 

We introduce a new field t) such that 

r(j){r, t) ’ 

where A is a constant [^. The field equations are 
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where overdots and primes denote derivatives with re¬ 
spect to t and r, respectively. The first three equations 
above correspond to the tr and 09 components of Eq. 
(^), and the wave equation, respectively. The final equa¬ 
tion is a convenient linear combination of the tt and rr 
components of Eq. (|^), and the wave equation. 

To find the time-independent solutions as van Putten 
does, all time derivatives appearing in Eqs. (i are set to 
zero, yielding 
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y ^v' = 2Z(^ + v'), 

r 

where Z = and is identical to that defined 

by van Putten ||p. Details regardin^he solution of this 
system of equations can be found in ; here we will only 
quote the results. The metric components are 

= 1 - AZ + -Z {1 - Z) = -{Zi + Z){Z2 - Z) 
e e 
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where e is an integration constant, and the constants Zi 
and Z 2 are given by 

2'i = 2e — i i \/l — 4e + Ifie^ 

^2 = -2e + i + i\/l-4e+ Ifie^. (9) 

The field Z is found from the transcendental equation 
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Note that the field ip can be recovered once Z, A, and v 
are known. 

As van Putten points out, Eq. ( 0 ) has four solutions, 
only one of which is Schwarzschild-like in its exterior (van 
Putten’s type la (^). For this solution we have e > 0 and 
Z —> Z 2 as r ^ 0 while Z ^ 0 as r ^ cx). 

It is worthwhile to consider the small r behavior of 
these fields. In terms of the integration constant e, this 
behavior is found to be 
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If in addition to small r, we consider the limit of small e, 
these expressions reduce to 
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Fig. ^ displays the solution to Eqs. (||) subject to the 
initial conditions derived from ( 0 ) for e = 1/100 (the 
same value shown in Q). 

Having constructed these static solutions, we can now 
address the question of their stability. For a given e, such 
a time-independent solution used in numerical relativity 
as an approximate black hole, were it not stable, would 
either collapse to a black hole or possibly disperse leaving 
flat space. 



z - ln(r) 

FIG. 1. Unperturbed fields for e = 0.01. As e ^ 0, the 
field Z approaches a step-function (see Fig. ^ and the field 
e'*' becomes more sharply peaked. This value of e is chosen to 
correspond with Fig. 2 of [1]. 


To investigate stability, we perform a standard linear 
perturbation analysis. As such, we consider the case that 
the fields do possess a small time dependent part and 
make the following expansion for small <5 

Ip —> -ipoir) + Sipi{r, t) 

V ^ vo{r) + 5vi{r,t) (13) 

A ^ Ao(r) 

where subscripts 0 and 1 denote the unperturbed and 
perturbed fields, respectively. 

We substitute the expansion (|^ into the full set of 
Eqs. ( 0 ) , keeping only terms to linear order. Because the 
unperturbed fields satisfy ( 0 ) by construction, we are left 
with the following linear equations for {tpi , vi, Ai) 
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We perform the standard Fourier decomposition of 
ipi{t,r) 


'0i(r, cr) = J e'''^*ipi{r,t)dt (15) 

and the other perturbed fields. On substitution into (0)’ 
the defining relation for -ipi then decouples from Ai and 
Di, yielding a single equation 
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It turns out that we need to use only the first couple of 
terms in the expansion to get accurate results. 
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which can be solved for the mode 4’i{r). We have 
thus reduced the perturbation problem to a single one¬ 
dimensional ODE with the undetermined eigenvalue tr^. 



eter e. The cross-marks and open circles represent data from 
two independent evaluations of the modes. Their correspon¬ 
dence indicates a quite small uncertainty. The vertical dotted 
line denotes the value of e for which Figs. &| and P are 
computed. 



FIG. 3. The unperturbed field Zo(r) for 0.025 < e < 0.24 
(uniform in t). The held Zo(r) approaches a step function as 
e is decreased. The approach of the unperturbed Zq to zero at 
large 2 ; signihes the approach to Schwarzschild in that region. 
For 2 < 0, however, the solution clearly is not Schwarzschild. 


Solution of Eq. requires appropriate boundary 

conditions on It is common to enforce regularity 

at the origin of a spherically symmetric spacetime, but in 
the current case the unperturbed solution is itself irregu¬ 
lar. Thus, assuming regularity of the perturbation might 
seem improper. However, although ^0 is logarithmically 
divergent at r = 0, exp{'ipQ) is regular at the origin, going 
to zero as a positive power of r. Hence, it is not un¬ 
reasonable to impose regularity on the field e'^. Further, 
because exp(^) = exp{ipQ+S'il^i) = exp('!/;o)(l-l-(5V’i), it is 
reasonable to assume the regularity of ipi at the origin. 
At most, the mode could have a logarithmic diver¬ 
gence, however, if we can find an unstable mode with the 
stricter criterion of regularity, then the solutions are still, 
in general, unstable. 

Enforcing the assumption of regularity of ipi, allows us 
to find a series expansion for ipi near the origin. For very 
small r, Eq. becomes 

0 ='ipi + + a'^CtpirP. (17) 

r 

The positive coefficient C is determined from Eq. ( 0 ) 
and depends on e. The exponent p = -^{Zi + Z 2 — 1) 
likewise depends only on e and in such a way that p > —2 
for e > 0. We can now find an expansion for r 1 
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Because of the linearity of the problem, '0i(O) can be 
arbitrarily chosen as it reflects the scaling in the problem. 


Given a background solution to Eqs. (^) for a particular 
e, we can now solve the eigenvalue problem (M) for the 
modes ipi (r) and corresponding characteristic frequencies 
In our particular case, the instability of the original 
soliton solutions is indicated by the existence of one or 
more exponentially growing modes. These are solutions 
to the perturbation equations with negative eigenvalues: 

< 0 . 

In practice, we integrate the unperturbed equations 
and the perturbation equation simultaneously from r « 0 
to large r, looking for a solution which has a negative 
eigenvalue and obeys the boundary conditions. Specifi¬ 
cally, we demand that the mode be finite at the origin 
and vanish asymptotically. We use a standard ODE in¬ 
tegrator and standard shooting techniques in our search. 

Although our search has not been exhaustive, we 
generically find precisely one growing mode for each value 
of e. This is sufficient for us to conclude that the static 
solutions are unstable. The eigenvalues found for these 
solutions are shown in Fig. |^. 

Having found the perturbation modes, looking at the 
limit e —> 0 is instructive. As this limit is approached, 
the unperturbed solution becomes more and more like 
Schwarzschild in the exterior, and this resemblance is 
precisely the reason why the family has been proposed 
as a good model of a black hole. With this in mind, 
one may wonder why Fig. ^ shows that, as e ^ 0, there 
is still a growing mode. Certainly these results do not 
show Schwarzschild to be unstable; rather we point out 
that within the “effective horizon” of this approximate 
black hole, the solution is very different from the interior 
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Schwarzschild solution for all e (see Fig. Hence, it is 
reasonable to assert that the solution is unstable for any 
e, including the solution e = 0. 
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FIG. 5. Comparison of modes computed in perturba¬ 
tion theory (solid lines) with modes computed from the full 
non-linear evolution (dashed lines). The non-linear modes 
were computed by taking a late-time profile from the evolu¬ 
tion and subtracting the initial data for that field, as defined 
in Eq. (p^). Perturbative modes have been numerically trans¬ 
formed to the Einstein frame. All helds are rescaled to the 
interval [0..1] and are plotted with respect to r. 


FIG. 4. Series of snapshots of the product rllj (related 
to the time derivative of ip) in the Einstein frame. A small 
perturbation {AM/M — 0.018%) at large radius is introduced 
to the initially static solution (e = 0.127). The perturbation 
passes through the singularity at r = 0 (between t = 24 and 
28) and escapes to r = oo. As the perturbation passes the 
red-shift horizon (as it propagates inwards, the perturbation 
is seen to experience a blue-shift), the excitation of a growing 
mode is clearly seen. 

Further, since it is the case that to an outside ob¬ 
server, the e = 0 solution is indistinguishable from that 
of Schwarzschild, it is logical to assume that any pertur¬ 
bation of the solution will not change the view of this ob¬ 
server. In other words, as e ^ 0, any perturbation should 
have decreasing support outside the “effective horizon” of 
the approximate black hole. Indeed, we observe this kind 
of behavior. As one decreases e, the profile of the mode 
ipi is seen to approach that of a delta function at the po¬ 
sition where the apparent horizon would asymptotically 
form. 

We find additional evidence for instability by evolv¬ 
ing the static solutions (with small perturbations added) 
using the full time-dependent equations of motion [^. 
This has the added benefit of providing information con¬ 
cerning the states to which the static solutions evolve 
when perturbed. To this end, we have adapted a pre¬ 
viously developed spherically symmetric code for Brans- 
Dicke theory which allows us to evolve these solutions. 
Because the evolutions in ||^ are performed in the Ein¬ 
stein conformal frame (as compared to the Brans-Dicke 
frame in which van Putten works and in which we have 
worked thus far), we transform the fields used above to 
the Einstein frame. In this frame, the field equations are 
equivalent to those for a massless scalar field minimally 
coupled to gravity. 


After recovering the static solutions in the Einstein 
frame and inputing a transformed solution into the code, 
we introduce a small ingoing perturbation to the fields 
at large radius. For generic values of e, we find that van 
Putten’s solution either collapses or disperses after the 
perturbation reaches the high-red shift horizon. Fig. ^ 
clearly demonstrates the instability for a specific e. In 
this case, the perturbation induces collapse to a black 
hole. 

In order to facilitate comparison of the perturbation 
results with those of the fully nonlinear evolution, we 
first define the quantities below in terms of the fields 
found in ||] 

Oi = a{r, t) — a{r, 0) 

$i = $^(r,t)-$^(r,0) (18) 

Bi = nj(r,t) -n5(r,0) 

such that fields with a bar and a subscripted 1 denote 
nonlinear deviations from the unperturbed solution {non¬ 
linear modes). By transforming the modes found from 
linear perturbation theory to the Einstein frame, we may 
now compare directly the linear modes with the nonlin¬ 
ear modes (Bi, oi, $i). In Fig. ||we show all the modes 
rescaled to the unit interval. From the near congruence, 
we conclude that we are observing the actual evolution 
of the growth of these perturbation modes. 

In addition to confirming our perturbative results, the 
full evolution provides evidence that these static solutions 
represent critical solutions to black hole formation. By 
this we mean that these solutions represent a boundary in 
the space of solutions between those that form black holes 
and those that do not. To demonstrate this criticality we 
begin with an unperturbed solution in the Einstein frame 
and add a perturbation to the fields. In this case instead 
of the arbitrary perturbation at large r shown in Fig. 
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we add the predicted mode found in the perturbation 
analysis (solid lines of Fig. ||). When this mode is added 
with some small positive amplitude, we invariably see 
collapse of the van Putten approximate black hole to a 
genuine black hole (see the solid lines in Fig. ||). In 
contrast, when the perturbation is added with a small 
negative amplitude, we see dispersal of the solution (see 
the dashed lines in Fig. Thus it would appear that 
this solution sits at the threshold between solutions that 
form black holes and those that disperse ||]. 



FIG. 6. Uniform-in-time series of solution when initially 
perturbed with the predicted mode {AM/M — 0.004%, 
e = 0.127). Shown here is the field II^ (related to the time 
derivative of tp). The evolution shown in solid results in col¬ 
lapse to a black hole. The dashed line shows the evolution 
resulting from switching the sign of the perturbation at the 
initial time. That the initial sign of the introduced perturba¬ 
tion separates eventual collapse from dispersal indicates that 
the unperturbed solution is critical. 


Having found evidence that these are threshold so¬ 
lutions, we are led to ask if they are attracting. For 
them to represent attracting critical solutions (interme¬ 
diate attractors) the unstable (relevant) mode which we 
find must be shown to be unique |^]. Because both the 
perturbation analysis and the full evolution appear to in¬ 
dicate the presence of only a single unstable mode (see 
Fig. |), we suspect that these solutions might well repre¬ 
sent an intermediate attractor for black hole formation. 
We plan to address this issue in future work. 

From the results presented here, these solutions would 
appear to be analagous to the n = 1 Bartnik-McKinnon 
(BM) solution in the Einstein-Yang-Mills (EYM) sys¬ 
tem 1^ . This static solution was found to be an interme¬ 
diate attractor in the gravitational collapse of spherically 
symmetric SU{2) fields with one side of the threshold be¬ 
ing black hole formation and the other dispersion of the 
Yang-Mills field [^j. 

After completion of this work, we became aware of 
other work which had considered solutions similar to 


those examined here. Static, spherically symmetric so¬ 
lutions to the minimally coupled Einstein-Klein-Gordon 
equations were studied by Buchdahl and later by 
Wyman [|m[. These solutions were written down in the 
Einstein frame in contrast to van Putten who works in 
vacuum Brans-Dicke (which is conformally equivalent to 
the Einstein-Klein-Gordon system). In addition, Jetzer 
and Scialom were able to show that Wyman’s solutions 
are generically unstable by establishing the existence of 
a negative upper bound for the lowest eigenvalue of the 
perturbation p^ . 
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